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Abstract 

A critical review of the theoretical understanding of the lifetimes of beauty hadrons 
is given. A model-independent analysis using the heavy-quark expansion allows 
for a description of the lifetime ratios t(B~)/t(B°) and t(Ai > )/t(B°) in a small 
region of parameter space. It is demonstrated that the lifetime ratio t(B~)/t(B°) 
cannot be used to extract the decay constant of the B meson. Implications for the 
semileptonic branching ratio and charm yield in B decays are pointed out. 
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1 Introduction 



The theoretical description of inclusive decays of hadrons containing a heavy quark ex- 
ploits two observations JJ: bound-state effects related to the initial state can be 
calculated using the heavy-quark expansion, and the fact that the final state consists of 
a sum over many hadronic channels eliminates the sensitivity to the properties of indi- 
vidual final-state hadrons. The second feature rests on the hypothesis of quark-hadron 
duality, i.e. the assumption that decay rates are calculable in QCD after a "smearing" 
procedure has been applied 0. In semileptonic decays the integration over the lepton 
spectrum provides a smearing over the invariant hadronic mass of the final state (so- 
called global duality), whereas for nonleptonic decays, where the total hadronic mass is 
fixed, the summation over many hadronic final states provides an averaging (so-called 
local duality). The calculation of the lifetimes of heavy hadrons relies on the stronger 
assumption of local duality. 

At present, quark-hadron duality cannot be derived from first principles, although 
some interesting attempts to address the problem have been discussed recently [§. 
The validity of global duality (at energies even lower than those relevant in B decays) 
has been tested experimentally using high-precision data on hadronic r decays [10]. The 
study of beauty lifetime ratios provides a test of local duality at energies of order a few 
GeV. Unlike the semileptonic branching ratio of B mesons, lifetime ratios are insensitive 
to possible new physics contributions and thus probe directly low-energy strong inter- 
actions. The common folklore is that the current experimental value r{h.^)/r{B) ~ 0.8 
cannot be accommodated by theory, indicating sizeable duality violations. If this were 
to be true, we could not trust the predictions of the heavy-quark expansion for other in- 
clusive quantities, such as the semileptonic branching ratio and charm yield in B decays, 
which are considered sensitive probes of new physics. A calibration of the theoretical 
uncertainties in the predictions for these quantities is therefore most important. 

Here we shall reconsider the theoretical description of beauty lifetime ratios taking a 
conservative point of view, which avoids making uncontrolled approximations We 
concentrate on the ratios r(S - )/r(S°) and r(A 6 )/r(S°), for which precise experimental 
data exist. Our model-independent approach is less ambitious than previous analyses in 
that it gives up some predictive power. Yet, we find it important to question commonly 
made assumptions before concluding a failure of the heavy-quark expansion. We find 
that at present it is still possible to reproduce the experimental data on beauty lifetimes 
in the framework of the heavy-quark expansion without invoking duality violations. Only 
a better control of theoretical uncertainties could help to clarify the situation. 

2 Heavy- Quark Expansion 

The inclusive decay width of a hadron Hj, containing a b quark can be written in the 
form 

T(H b - X) = Im (H b \ T \H b ) , (1) 

M H b 
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where the transition operator T is given by 

T = *|d 4 xT{£ eff (x),£ cfr (0)}. (2) 



Here £ e ff is the effective weak Lagrangian renormalized at the scale QT2J, [13|| . The 
leading contributions to the transition operator are given by the two-loop diagrams shown 
on the left-hand side in Fig. |l| Because of the large mass of the b quark, the momenta 
flowing through the internal propagator lines are large. It is thus possible to construct an 
operator product expansion (OPE) for the transition operator, in which it is represented 
as a series of local operators containing the 6-quark fields. The operator with the lowest 
dimension is bb. It arises by contracting the internal lines of the first diagram. The only 
gauge-invariant operator with dimension 4 is bipb; however, the equations of motion 
imply that this operator can be replaced by mj)b. The first operator that is different 
from bb has dimension 5 and contains the gluon field. It arises from diagrams in which 
a gluon is emitted from one of the internal lines, such as the second diagram shown in 
Fig. [I]. From dimension 6 on, a large number of operators appears. For dimensional 
reasons, the matrix elements of higher-dimensional operators are suppressed by inverse 
powers of the 6-quark mass. Thus, the total inclusive decay rate (i.e. the inverse lifetime) 
of a hadron Hf, can be written as @. U 

G 2 F ml\V cb \ 2 \ - (bg s a, v G^b) {n) (0< w) ) \ 

where the prefactor arises from the loop integrations, q are calculable coefficient func- 
tions, and (O) are the (normalized) forward matrix elements between H b states. These 
matrix elements can be systematically expanded in powers of 1/ m^ using the heavy-quark 
effective theory (HQET) 0. The result is @, |] 

m = l _ Aw-^c.m +0(l/mlh 

(bg.a^Cb) = 2tf-j(H b ) + 0(l/m„) , (4) 

where ^(H^) and n 2 G (Hb) are the matrix elements of the heavy-quark kinetic energy 
and chromomagnetic interaction inside the hadron Hf,, respectively | 15| . Inserting these 
relations into (^) and taking the ratio of two lifetimes yields |JTTJ 

t(B-) 



- -l + 0(l/m bJ , 



3\ 



= (0.98 ±0.01)+ 0(l/ml), (5) 



where cq = —\ — Ic^jcj, ~ 1.2, and isospin symmetry is assumed for the matrix ele- 
ments of B~ and B° mesons. The hadronic parameters appearing at order l/m| in the 
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Figure 1: Leading contributions to the transition operator T (left), and the correspond- 
ing operators in the OPE (right). The open squares represent a four-fermion interaction 
of the effective weak Lagrangian. 



ratio t(Ai)/t(B°) can all be determined from spectroscopy. The matrix elements of the 
chromomagnetic interaction are given by 



/4(£) = -(Ml* - M%) « 0.36 GeV 2 







(6) 



where the second relation follows from the fact that in the ground-state A& baryon the 
light constituents have total spin zero. The individual matrix elements of the kinetic 
energy are not physical quantities because of renormalon ambiguities [IB"]; however, their 
difference obeys the relation 



/4(A*) - A4(2J) 



M B M D ( M\. — Ma 



Mb -Af; 



D 



3 M B * - M D 

4 M B - M 



D 



(7) 



which implies /i 2 (A b ) = fJ%(B) up to small corrections of order 0.01 GeV 2 . 

The mo del- independent theoretical predictions in (M) may be compared with the 



average experimental values for the lifetime ratios, which are [17| 



r(B-] 
t(B°) 



1.06 ±0.04 



r(A b ) 
t(B°) 



0.79 ±0.05. 



(8) 



Whereas the first ratio is in good agreement with the theoretical expectation, the low 
value of the lifetime of the baryon constitutes a puzzle. To understand the fine struc- 
ture of the lifetime ratios requires to go further in the l/m& expansion. At order l/m\ 
in the expansion of nonleptonic decay rates there appear four-quark operators whose 
matrix elements explicitly depend on the flavour of the spectator quark(s) in the hadron 
Hb, and hence differentiate between hadrons with different light constituents. These 
nonspectator effects receive a phase-space enhancement factor of 0(16tt 2 ) with respect 
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to the leading terms in the OPE, which makes them more important than the corrections 
of order l/m%. This can be seen from Fig. ^, which shows that the corresponding con- 
tributions to the transition operator arise from one-loop rather than two-loop diagrams. 
In total, a set of four four-quark operators is induced by nonspectator effects. They are 

O x = 67^(1 -75)9^(1 -75)6, 

2 = 6(1-75)99(1 + 75)6, 

Tx = 67„(1 - 75Kgg7 M (l - 75)46, 

T 2 = b(l- l5 )t a qq(l + l5 )t a b } (9) 

where q is a light quark, and t a are the generators of colour SU(3). Since hadronic 
matrix elements of four-quark operators are notoriously difficult to calculate, one needs 
to parametrize them. For the S-meson matrix elements renormalized at the scale rrif,, 



we define 11 



(B\ \B), =mb = B t f l B M B , (BIT, \B), =mb = e t f B M B , (10) 

where /# is the decay constant of the B meson. The values of the dimensionless hadronic 
parameters Bi and £j are currently not known. Ultimately, they may be calculated using 
some field-theoretic approach such as lattice gauge theory. Some (but not all) of these 
parameters may also be extracted from a precise measurement of the lepton spectrum 
in the endpoint region of semileptonic B decays fll8|. We note that in the vacuum 



saturation (or factorization) approximation |I9| Bi = 1 and = at some scale fx, 
where the approximation is assumed to hold. More generally, the large-iV c counting 
rules of QCD imply 

Bi = 0(l), ei = 0{l/N e ). (11) 

The theoretical expressions can also be written in terms of parameters Bi(fi) and 
Ei(fjj) renormalized at a scale \x 7^ m&, using evolution equations for the hadronic matrix 
elements in the HQET. Such a rewriting does not induce any uncertainty, as it is simply 
a reparametrization of the results. Still, it may be appropriate to use a lower renor- 
malization point if model calculations of the parameters Bi and £j are used. To give an 
example, we relate the parameters in (|i~0|) to ones renormalized at a low hadronic scale 



A*had chosen such that a s (/ihad) = 0.5. In defining these parameters, we also renormalize 



the decay constant /# at the scale //had- Then the relations are |TT 

Bi w l.OLBi(^had) - 0.24£i(/i had ) , 

e t w -0.05Bi(fu had ) + 0.73£i(/i had ) . (12) 

Next consider the matrix elements of the four-quark operators between A& baryons. 
Heavy-quark spin symmetry implies the relations (O2) = — \{Ox) and (T 2 ) = — \(Tx) [PH , 
which reflect the fact that the light degrees of freedom in a baryon are in a state with 
total spin zero. This leaves two independent matrix elements of the operators Ox and 
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b Tj qqTjb 





Figure 2: Nonspectator contributions to the transition operator T (left), and the 
corresponding operators in the OPE (right). I\ denote combinations of Dirac and colour 
matrices. 



7 5 ) b l instead 
-10, + \O x . 



T\. It is convenient to introduce the operator 0\ = &*7 M (1 — 75)^ ^7^(1 
of Ti are colour indices) by making use of the Fierz identity 7\ = 
The analogue of the vacuum saturation approximation in the case of baryons is the 
valence-quark assumption, in which the colour of the quark fields in the operators is 
identified with the colour of the quarks inside the baryon. Since the colour wave function 
of a baryon is totally antisymmetric, the matrix elements of 0\ and 0\ differ in this 
approximation only by a sign. Hence, it is natural to define 



(A b \Oi\A b ) 



ti=m b 



-B (A b |Oi|A 6 ) 



fi=m b 1 



(13) 



with B = 1 in the valence-quark approximation. For the baryon matrix element of 0\ 
itself, we adopt a parametrization guided by the quark model and write 



1 



2 Ma, 



(A 6 | Ot \A b ) 



B 



fj,=m b 



48 



(14) 



In the quark model r is the ratio of the squares of the wave functions determining the 
probability to find a light quark at the location of the b quark inside the A^ baryon 
and the B meson |]20], 21| . The parameter r is by far the most uncertain one entering 
the theoretical formulae for beauty lifetimes. Theoretical estimates range from r ~ 0.1- 
0.5 [E2l |23| to r pa 1-2 p4[. In the quark model, r can be related to the ratio of the spin 



(15) 



splittings between E& and baryons and B and B* mesons, leading to 



4 Ml, 



Ml 



3 M|» 



Ml 



Using the (still preliminary) result Ms* —M^ b = (56 ±16) MeV reported by the DELPHI 
Collaboration ||25|| , one then obtains r ps 1.8 ± 0.5. 
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3 Model-Independent Analysis of Lifetime Ratios 



Including the nonspectator contributions of order 1/mf allows us to refine the results in 
(H) at the prize of introducing several unknown hadronic parameters. We obtain [] I 



t(B- 



t(B°) 

t(A») 
t(B°) 



1 + £ k i B 1 + k 2 B 2 + k 3 e 1 + k±e 2 



0.98 + £ p x B t + p 2 B 2 + p 3 e 1 + p 4 e 2 + (p 5 + p 6 B)r , (16) 



where £ = (/b/200 MeV) 2 . The coefficients fcj and pt depend on the quark-mass ratio 
m c /m h and the values of the Wilson coefficients c±{m h ) appearing in the effective weak 
Lagrangian. They also depend on the perturbative scheme adopted in the calculation. 
So far, we have discussed the scheme in which the standard QCD evolution is used 
to obtain the effective weak Lagrangian at the scale mj, where the matrix elements in 
the heavy-quark expansion are evaluated. To investigate the perturbative uncertainties, 
one may also renormalize the weak Lag rcLIlgicLll 8bt cL SCctl.6 fi = Krrib with k = 0(1), then 
perform the heavy-quark expansion, and finally evolve the hadronic matrix elements from 
nrrib to 777.6 using evolution equations in the HQET. If it was not for the truncation of 
perturbation theory, the two procedures would yield to the same result. The dependence 
of hi and pi on the choice of k is therefore a measure of the truncation error. In Table [3], 
we give the values of these coefficients for the three choices k = 1/2, 1 and 2. Because 
of an accidental cancellation between the contribution from different operators in the 
effective weak Lagrangian, the coefficients k\ and have a significant dependence on 
the choice of the "matching" scale, whereas the other coefficients are either negligibly 
small or stable with respect to variations of k. 



Table 1: Coefficients fcj and pi appearing in the predictions for the lifetime ratios. We 
use c + (m&) = 0.86, c^(rrib) = 1.35, and m c /mb = 0.29 



K 


h 


k 2 


h 


fc 4 






1/2 


+0.046 


0.003 


-0.737 


0.201 






1 


+0.021 


0.004 


-0.699 


0.195 






2 


-0.007 


0.007 


-0.666 


0.189 








Pi 


P2 


P3 


Pa 


P5 


P6 


1/2 


-0.003 


0.003 


-0.178 


0.201 


-0.017 


-0.023 


1 


-0.003 


0.004 


-0.173 


0.195 


-0.012 


-0.021 


2 


-0.006 


0.007 


-0.168 


0.189 


-0.008 


-0.020 



Previous analyses of beauty lifetimes [^3, |27| have avoided the proliferation of un- 
known hadronic parameters by making strong model assumptions. The meson matrix 
elements of four-quark operators have been estimated using the vacuum saturation ap- 
proximation at a low hadronic scale /ihad which, according to (0), yields Bi w 1.0 and 
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e i pa —0.05 at the scale m&. The baryon matrix elements have been estimated adopting 
the valence-quark approximation B = 1 and assuming r = 0(1). These assumptions 
lead to the model predictions 

= 1 + (0.05 ± 0.03) x ( , 4^0.95. (17) 
t(B°) K ' V200MeV/ r(B°) K J 

Whereas the first relation is in good agreement with the experimental result in (§), the 
predicted value for the ratio r(Ab)/r(5°) is much larger than the experimental one. This 
discrepancy is referred to as the "A& lifetime puzzle" . 

It has been advocated in Refs. p7| that a precise measurement of the lifetime ratio 



t(B )/t(B°) would allow for a determination of the 5-meson decay constant /g — a 
proposal which has been taken seriously by the DELPHI Collaboration P5[. Note that 



the ambiguity in the choice of the matching scale alone implies a large uncertainty in the 
coefficient in front of j\ in QT7D, which is ignored in Refs. [J27| . In addition, the numerical 
results for the coefficients ki in Table [I] show that assuming exact vacuum saturation is 
not justifiable in the case of the ratio t(B~)/t(B°). To neglect in (|16|) the deviation 
As\ of the unknown parameter e\ from the value —0.05 predicted by factorization at a 
low scale would only be justified if |Aei| <C 0.1. The large- N c counting rules of QCD do 
not allow us to impose such a strong restriction. 

Our goal here is to investigate in a model-independent way how the "predictions" in 
(|T7|) are altered if model assumptions are avoided. To this end, we allow all hadronic 
parameters to vary within reasonable limits and use the large- Af c counting rules and 
other theoretical prejudices only to restrict the allowed parameter ranges. To be spe- 
cific, we take random distributions inside the following ranges: Bi, B G [2/3,4/3], 
£i G [—1/3, 1/3], r G [0.25, 2.5], and £ G [0.8, 1.2]. For each parameter set, the lifetime ra- 
tios t(B~)/t(B°) and r(A b )/r(5°) are computed for the two perturbative schemes with 
k — 1/2 and 1. The results are shown in Fig. || In general, the scheme with k — 1/2 
gives results closer to the experimental data. There are two important observations: 
first, only a small tail of the theoretical distribution for the ratio t(A&)/t(.B ) reaches 
into the region preferred by the experimental data; secondly, the predictions for the ra- 
tio t(B~)/t(B°) have an almost flat distribution between 0.85 and 1.25, indicating that 
without a better knowledge of the hadronic parameters it is not possible to predict this 
ratio with good precision. At present, not even the sign of the deviation from unity can 
be predicted without model assumptions, in stark contrast with the model-dependent 
result in fllTD. 

In the next step, we focus on the parameter sets which give results within the 2a 
ellipse around the central experimental values. For these, the projected distributions 
for all hadronic parameters are shown in Fig. |j. The distributions of B\ and B 2 and, 
to a lesser extend, B and £ are rather flat, indicating that these parameters are little 
restricted by the data. The remaining three parameters are highly constrained, however. 
The parameter ei, in particular, is determined by the experimental data within a rather 
limited range: S\ ~ — (0.1 ± 0.1), in good agreement with the expectation E\ ~ —0.05 
based on the vacuum saturation approximation. The distributions of the parameters e 2 
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Figure 3: Theoretical results for the lifetime ratios for random choices of the hadronic 
parameters within the limits specified in the text. The ellipses in the upper plots show 
the la and 2a contours around the central experimental values. In the lower plots, the 
projections on the individual lifetime ratios are shown. Full lines correspond to k = 1, 
dashed ones to k = 1/2. 



and r are strongly correlated. To reproduce the experimental data requires to have 82 
negative and r > 1. The correlation between the two parameters in shown in the last plot 
(for the scheme with k — 1/2). If we assume that \e2\ < 0.15, for instance, then r must be 
larger than about 1.25, which is higher than most expectations [22], |23[. To summarize, 
although some large hadronic parameters are needed to explain the observed pattern of 
beauty lifetimes, a conventional explanation of the A& lifetime puzzle cannot be excluded 
before reliable field-theoretic calculations of the matrix elements of four-quark operators 
become available. 



4 Conclusions and Implications 

The heavy-quark expansion, supplemented by the assumption of quark-hadron duality, 
provides the theoretical framework for a systematic calculation of the lifetimes and in- 
clusive decay rates of hadrons containing a heavy quark. Whereas this formalism works 
well for the description of inclusive semileptonic decays, two potential problems related 
to nonleptonic decays have become apparent in recent years: the low experimental value 
of the lifetime ratio r(A^)/r(i? ), and the low values of the semileptonic branching ratio 
and charm yield in B decays. 

In order to obtain a detailed understanding of beauty lifetimes it is necessary to go 
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Figure 4: Distributions of the hadronic parameters for the simulations giving 
within 2a of the central experimental values. 
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to order l/m| in the heavy-quark expansion, at which the matrix elements of four-quark 
operators appear. They describe the physics of nonspectator effects, i.e. contributions 
in which a light quark in a beauty hadron is actively involved in the weak interac- 
tion. We have presented a model-independent study of such contributions introducing 
a minimal set of hadronic parameters, which eventually may be determined using some 
field-theoretic approach such as lattice gauge theory. We find that in I?-meson decays 
the coefficients of the colour-octet nonfactorizable operators are much larger than those 
of the colour-singlet factorizable operators, and therefore the contributions from the non- 
factorizable operators cannot be neglected. We also find that the theoretical predictions 
for lifetime ratios still suffer from large perturbative uncertainties, which could only be 
reduced if the short- distance coefficients in the heavy-quark expansion were calculated 
to next-to-leading order. The theoretical prediction for the ratio t(B~)/t(B°) is in 
agreement with experiment; however, our present ignorance about the matrix elements 
of four-quark operators does not allow us to calculate this ratio with an accuracy of 
better than about 20%. As a consequence, it is not possible to extract the 5-meson 
decay constant fs from a measurement of t(B~)/t(B°). 

The short A b lifetime remains a potential problem for the heavy-quark theory. If the 
current experimental value persists, either some hadronic matrix elements in the baryon 
and/or meson sector must be larger than naive expectations, or local quark-hadron 
duality fails in nonleptonic inclusive decays. In the latter case, the explanation of the 
puzzle of the A& lifetime lies beyond our present capabilities. We stress, however, that 
in view of our results it is premature to conclude a failure of the heavy-quark expansion 
caused by sizeable duality violations. There is a small region of parameter space where 
the experimental data can be accommodated by theory. 

If indeed the observed pattern of beauty lifetimes is caused by a conspiracy of some 
large hadronic matrix elements, which are the implications and tests of this scenario? 
A crucial test will be to prove (or disprove) with a field-theoretical calculation that at 
least one of the two parameters £2 and r has the large value required to explain the 
data. At present, efforts to compute these parameters using lattice gauge theory are 
under way |[29|| . An important implication concerns the theoretical understanding of the 
semileptonic branching ratio and charm yield in B decays. These issues are discussed at 
length elsewhere in these Proceedings pOf . All theoretical calculations of these quantities 



so far have neglected nonspectator contributions [ 31 1 — [ 33 1 . The most recent theoretical 



predictions are [[LI 



f 12.0 ±1.0%; /i = m b , 
SL ~ 1 10.9 ± 1.0%; ft = m b /2, 

f 1.20 =f 0.06; = m b , 
Hc ~ \ 1.21 =f 0.06; fi = m b /2, 

with correlated errors. Notice that the semileptonic branching ratio has a stronger scale 
dependence than n c . By choosing a low renormalization scale, values Bsl < 11% can be 
accommodated. Currently, the experimental values of Bsl and n c are still controversial; 
however, the value of n c tends to lie below the theoretical predictions. 
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Nonspectator contributions can change the predictions for the semileptonic branching 
ratio significantly; however, they have a negligible effect on n c . Irrespective of the precise 
values of the hadronic parameters, there is a model-independent relation An c « A5sl ~ 



0-1% pTT 1 , i.e. the relative effect on B$l is ten times larger than that on n c . A large 
negative value of 62 implies a positive contribution to -Bsl, which is not preferred by the 
data. This restriction favours the solution where the A& lifetime puzzle is solved by a 
large value of the baryon parameter r. 

The fact that the charm yield n c is very little affected by nonspectator effects is in- 
teresting. To reduce n c significantly below the values quoted in fll8|) requires an anoma- 
lously large branching ratio for charmless nonleptonic B decays. At present, a value 
of B(B — > X s + no charm) several times larger than predicted by the Standard Model 
is not excluded by the data. Such interesting scenarios are explored elsewhere in these 



Proceedings [34 
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